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FOREWARD 


This report presents the results of work performed by Mr. 
Jerry J. Rubli at Louisiana State University under contract 
NAS9-10464. This is the final report for the work done on the 
separation of the booster and orbiter under the Graduate Engine- 
ering-Practice in Mechanical Engineering program and it covers 
the analysis performed during 1970. The faculty advisor was 
Dr. Mehdy Sabbaghian and the program director was Dr. L. R. 
Daniel, Jr. Associate Directors were Mr. Charles Teixeira 
(NASA-MSC) and Dr. R. W. Courter. The principal advisor was 
Dr. P. H. Miller. 
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SUMMARY 


The analysis of the weight for four separational methods 
is presented. The computer program listing and results, which 
prove that it is feasible to design a pyrotechnic ram are 
enclosed. Finally, the energy versus weight curves for the 
four different techniques studied are presented. 
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number of spring coils 

universal gas constant for a specific gas, ft-lbf/ 
lbm°R 

1545 ft-lbf /lb mole °R 
mass fraction 

maximum yield strength for a helical compression spring, 
Ibf 

temperature, °R 
time, sec. 

internal energy, BTU/lbm 
gas velocity, in/s 
volume, in^ 
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specific volume, in /lbm 
weight, Ibf 

relative separational displacement, ft 
relative separational velocity, ft/s 

. 2 ' 

relative separational acceleration, ft/s 
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specific weight, Ibf /in 

spring deflection, ft 

booster rotation displacement, rad. 

booster rotational velocity, rad/s 
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booster rotational acceleration, rad/s 
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density, lbm/in 
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CHAPTER X 


INTRODUCTION 

Stage separation for space vehicles usually involves stages 
whose longitudinal axes are coaxial. The space shuttle envi- 
sioned by the NASA-Manned Spacecraft Center (MSC) in Houston, 
Texas, will involve two stages whose longitudinal axes will be 
parallel, but not coaxial. 

The space shuttle is a reuseable, two stage, manned space 
ship. Some of the possible missions of the space shuttle in- 
clude carrying men and supplies to an orbiting space station, 
satellite maintenance, and space rescue. The first stage of 
the shuttle, called the booster, will have a mass of 321,000 
lbm at separation and will be 203 feet in length. The second 
stage of the shuttle, called the orbiter, will have a mass of 
450,000 lbm at separation and will have an overall length of 
123 feet. Both vehicles will have the ability to land at a 
conventional airport following mission completion. The orbiter 
is attached to the booster in a "piggyback” arrangement such 
that the center of gravity of the orbiter is forward of the 
center of gravity of the booster. The shuttle is launched in 
a verticle attitude with the booster engines supplying the 
thrust. The booster stage engines are operated at full thrust 
level from lift off until a longitudinal acceleration of 2.5 g's 
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is obtained. The engines are then throttled to maintain the 
acceleration at a 2.5 g level. 

Planned booster/orbiter separation occurs at 267,000 feet, 

at which time the shuttle will have obtained a velocity of Mach 

2 

11 and the dynamic pressure will be 1 lb /ft . Emergen cy abort 
before the shuttle has attained an altitude of 130,000 feet is 
nearly impossible. From 130,000 feet to 230,000 an abort at 
separation followed by down range land recovery is possible. 

At and above 230,000 feet abort to orbit is possible. 

Since the overall weight of the shuttle influences the 
magnitude of the payload weight and the operational feasibility 
of the entire shuttle concept, the weight required for each 
separation technique is regarded as the main criterion in the 
selection of the separation method to be employed at staging. 

As the translational velocity due to separation has not 
yet been specified, the energy required for separation is not 
known. Therefore, a logical approach to be followed involves 
determining the equation which relates hardware weight and its 
energy delivery capability. One difficulty which does arise, 
however, is the fact that this functional relationship can some 
times be expressed only implicitly. 

The purpose of this analysis is to graphically express the 
energy capability versus hardware weight for the four separa- 
tion techniques. 
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These techniques include: 

1. separation rockets 

2. spring loaded ram 

3. non-ignitable gas ejectors 

4. pyrotechnic ram 



CHAPTER II 


SEPARATION ROCKETS 

The assumptions made are the following: 

1. 2 rockets are located on the booster only, and the 
resultant force acts through its center of gravity. 

2. The rocket thrust is instantaneously developed, 
remains constant, and is instantaneously stopped. 

3. "State of the art" upper and lower values of the vacuum 

specific impulse for the propellant, I , are 300 and 

sp 

275 lbf-sec/lbm 5 , respectively, (according to NASA 
personnel) . 

4. The range of values for the mass fraction R m £> is from 
.85 to .90. The mass fraction is defined as the mass 
of the propellant divided by the total rocket mass. 

Referring to Figure (II- 1), the following dynamical relations 
govern the relative motion between the booster and orbiter during 
the operation interval. 
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beg 


fb 

“b Sc 


(II- 1) 
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x = — F.dt + c. 

beg J h 1 

Vg = J [$I F b dt + c l] dt + °2 


where 


and C 2 are integration constants 

Since the force F, is constant with time we can integrate 

b 

Equations (II-l), (II-2), (II-3) from the initiation of separa- 
tion until separation has been completed in order to obtain the 
final acceleration, velocity, and displacement 


• • 

X, 


^ g c 

% c 


(II-4) 


X„ = 


Vf 

■"b 
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X. = 
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2m b s c 


(II-6) 


The energy delivered to the booster/orbiter system is given 


as 


(n-7) 


E = 


* 

l 



dx + c 


where 


c, is an integration constant. 

Since is not a function of displacement, one can evaluate 
Equation (II-7) from the beginning of separation until the final 
displacement is obtained This provides 


E = F b X f 


(II-8) 


Substituting from Equation (II-6) into (II-8) yields 



(II-9) 


After rearrangement the above equation becomes 


F b t f 



(11-10) 


The total impulse on the booster is defined as the product 
of the thrust acting on the booster and the time during which it 
acts on the booster. This may be expressed as 
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tot 


F t 
b f 


( 11 - 11 ) 


Since the product of the specific impulse and the rocket 
fuel mass is the impulse capability of one rocket engine, the total 

i 

impulse capability of both rocket engines is 


I =2 I m _ 
tot sp fuel 


(11-12) 


Using Equations (11-10) and (11-11) the total impulse may 
be expressed as 


tot 


2a b E 

s„ 


(11-13) 


nu , 
fuel 


Substituting this into Equation (11-12) and solving for 
gives 



(11-14) 


From the definition of R . as defined in assumption 4 at the 

mt 

beginning of this chapter, it can be shown that 


m 


rocket 


m_ . 
fuel 

R faf 


(11-15) 


Substituting into Equation (11-15) the expression for 
from Equation (11-14) yields 



m rocket 21 R 

sp mf 


(11-16) 


Equation (11-16) indicates that the total rocket mass for a 
constant thrust engine is a function of the mass of the booster, 
the energy imparted to the booster, the specific impulse of the 
fuel, and the mass fraction. If it is desired to obtain the 
mass of tandem engines simply multiply by two to obtain 



m 2rockets I R c 

sp mf 


(11-17) 


In the English Engineering system of units 


(Force in lbf) = (mass in lbm) 

8 C 


(11-18) 



Therefore, 


JH 

2rockets 2rockets g 


(11-19) 


Substituting for H^ro^ets ^ roin Equation (11-17) into Equation 
(11-19) leaves 


W 


2rockets 



I R r 
sp mf 


£L 




( 11 - 20 ) 


By using the lower values for the mass fraction, R mf , and 

the specific impulse, I , one can plot weight versus energy to 

sp 

obtain a lower bound for heavier and relatively inefficient 
rockets for the rocket engine design. Employing maximum values 
of these parameters will provide the upper bound for the lighter 
and more efficient engines. These curves are plotted in Figures 
(VI-1) and (VI-2). 


CHAPTER III 


SPRING LOADED RAM 

The location of the spring loaded ram would be directly below 
the c.g. of the orbiter but forward of the booster c.g. The 
approach taken in the analysis of the required weight for this 
technique is to determine the relationship between weight and 
energy for the helical compression spring only. Since it will be 
demonstrated that the weight of the spring alone is excessive, the 
spring loaded ram hardware need not be considered. 

The assumptions made in the analysis of the spring loaded ram 
technique are as follows: 

1. A single spring device is located on the booster structure 

2. The spring is a helical spring. This assumption was made 

r 

because for a helical spring the spring constant is 
constant. This will result in a much smaller initial 
force relative to final force than would be the case for 
a variable spring constant. NASA personnel have advised 
that the initial force relative to final force should be 
as small as is possible. 

According to Reference (1), the deflection of a helical 
compression spring is given by 
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The elastic energy stored in a spring is 


E = 


F6 


(III- 2) 


Substituting for 6 from Equation (III-l) and solving for F 
results in 


F = 


2EGd 12 


8c 3 n 


(III- 3) 


The maximum stress induced in a spring of this type is given 


by 


8KFc 

’s ” a 1 
nd 


(III-4) 


where 


S = maximum shear stress lbf/in 
s 


K = Wahl factor 


= 0 


4c - 1 0.615 " 

4c - 4 c 


D 


c = 


Replacing F by the expression given in Equation (III- 3) yields 


8K 


S = 
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TTd 


\ [2EGd 12 
18cn _ 4K 
tt 


EG 


12 


cd n 


(III-5) 


or 
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[cd 3 n ] = m GE 


y-y 


(III- 6) 


The weight of a helical spring is 


W . = Y (nD)n 

spring ' 4 


(III-7) 


Substituting for D - cd in the above equation results in 


Yrr r 3 I 

W , = — r~ I cd n 

spring 4 L J 


(7TI-8) 


From Equations (III-6) and (III- 8) one would obtain 


48yG 


W 


spring 


K 2 E 


(III-9) 


Since it is desired to get the lightest spring possible, one needs 

2 

to select the minimum values for K and yG/ S g . As c becomes large, 

K approaches 1 asympototically. Therefore, K * 1 is used as an 

upper limit since no spring could be designed which would be 

lighter than one corresponding to this condition. 

2 

In order to minimize ( yG/S ), one may notice that for most 

s 

a 6 2 

spring steels, y ™ *285 lbf/in . Also, G s 10.5 x 10 lbf/in for 

6 2 

hot wound springs and G = 11.5 x 10 lbf/in for most cold drawn 
springs. The maximum allowable value for S g is shown in Figure 


(XXI- 1) to be equal to 130,000 psi at d *> .5 in. The spring 
material is hot wound SAE G150 or SAE 9260. Using these values and 
a safety factor of 1.5 in Equation (III- 9) yields 


W 

spring 


<• 


48( .285 (l0.5 x 10 6 •i~)(l) 2 (1.5) 


. 3/ 
in 


'E 


in 


(130, 000) 2 


in 


(III- 10) 


Performing the indicated operations in Equation (III- 10) leaves 


W . = .0191 E (III- 11) 

spring 

Equation (III- 10) is the desired relation between energy 
capability and the design weight. As can be seen in Figure (VI- 1), 
the excessive weight of the spring alone rules out the necessity 
for designing the hardware required to retain the spring. 



9 1.0 1.1 1.2 1.3 1.4 1 , 

Bar Diame.ter (inches) 


FIGURE III-l 




CHAPTER IV 


GAS EJECTOR 

This device utilizes the thrust provided by the high speed 
ejection of a non-ignitable .gas* Therefore, there is no 
impingement of high temperature gases on the arbiter. Because a 
one pound increase in the weight of the orb iter results in an 
overall shuttle weight increase of six pounds, it is more feasible 
to locate this device on the booster. Since a minimum of booster 
angular motion is desired at separation it is practical to locate 
the separational thrust through the booster c.g. If two ejectors 
are located on the booster wings, one on either side of the fuselage, 
the booster will experience no rotation during separation. 

In analyzing this separation technique the following assumptions 
are made : 

1. The thrust, stagnation pressure, and stagnation temperature 
are constant. Since it will be later demonstrated that 
even a constant pressure results in excessive weight, the 
actual device would be unfeasible. 

2. The hydrogen or oxygen gas contained in the booster fuel 
tanks at normal tank pressure can be used as the high speed 
ejector gas, depending upon which one is more feasible. This 
assumption eliminates the weight associated with having 
separate tanks for the ejector. 


16 



17 


3. An isentropic process occurs across the ejector exit 
nozzle. This assumption allows the use of Equation (IV- 1) 
for the theoretical specific impulse. 

4. The gases exit at the ambient pressure at separation altitude, 
i.e., complete- expansion. This assumption was also necessary 
in order that Equation (IV- 1) could be utilized. 

These assumptions all result in a lighter weight for the gas 
being ejected. Since it will be later demonstrated that this 
technique results in an excessive weight, use of these 
simplifying assumptions is valid. 

According to Reference (2), the theoretical specific impulse 
of a gas flowing through a nozzle is given as 


sp 



7 * (■ - © 

m 0 


k-1 
k 


(IV- 1) 


where 

Pq is the stagnation pressure 
is the ambient pressure 
T q is the stagnation temperature 

Using this theoretical value will result in a more conservative 

(lighter) value for the gas weight. 

Substituting for R and g and evaluating gives 

c 


' sp =< 6 - 94 > \ 2 <irr-i> f ( l - 


(IV-2) 


:) 
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Since the largest value for the terms in the last parentheses 
is 1.0 for = 0, a higher value for the theoretical specific 
impulse, I* , will result if is set equal to zero. This will 
have the effect of decreasing the weight of required gas. Since 
it will be later shown that the gas ejector concept results in an 
excessive weight even with this assumption, it is valid to consider 

p 2 = 

Since the theoretical specific impulse is defined as the 
impulse divided by the mass of gas producing the thrust, we can 
obtain the required mass of gas for one ejector by dividing the 
impulse of one engine by the theoretical specific impulse. Since 
the impulse provided by one of the two gas ejectors is one half 
of the total impulse aeti-ng on the boos ter j 


m 

gas 



(IV- 3) 


From assumption 1, since the thrust is considered constant. 
Equations (II- 10) and (II- 11) give for this condition 



(IV-4) 


Substituting I tot from Equation (IV-4) into Equation (IV- 3) 


leaves 


19 


2m b E 

8„ 


m 


gas 
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(IV- 5) 


sp 


The mass of gas required for both ejectors is twice that as 
expressed in Equation (IV-5). It is given as 

X E 


8 . 


m 


2 gas 


(IV- 6) 


sp 


Substituting for I* from Equation (IV-2) into Equation (IV-6) 

sp 

leaves 


m 0 = .144 
2 gas 




o s c 


E 


(IV- 7) 


In deciding whether oxygen or hydrogen should be used as the 
ejecting gas, the gas which gives the minimum product of 
will result in less gas weight required, as can be seen in 


Equation (IV-7). For oxygen this product is equal to 4.42 x 10 


-2 


lbm 


lb mole °R 


. For hydrogen this product is equal to 1.39 x 10 


-2 


lbm 


lb mole °R 
gases . 


. Therefore, hydrogen is the more feasible of the two 


If in Equation (IV-5) this value of 1.39 x 10“^ ..w. 

- lb mole ^RL 

is substituted for y 0 > the mass of the booster , 

which is equal to 321,000 lbm, is substituted for m^, and a safety 
factor of 1.5 is applied, the mass of gas required to supply two 
high speed non-ignitable gas ejectors becomes 


= 2.53 \|e 


(IV- 8) 


m 


2 gas 


Applying Equation (11-18) to Equation (IV-18) , the weight 
for two ejectors can be expressed as 


W 


2 gas 


= 2.53 


v 


E 


J£_ 


(iv- 9) 


This is the desired relationship ^between energy and weight. This 
equation is graphically displayed in Figure (VI-1). Since the 
gas weight for this device is excessive, the- other components of 
such a' system need not be considered. 



CHAPTER V 


PYROTECHNIC RAM 

This chapter is divided into five sections. The first sec- 
tion deals with the description of this device and its operation. 
The second section is an analysis of the dynamics and energy in- 
volved in separating the two vehicles. The third section is con- 
cerned with the thermodynamics of the ram operation. The fourth 
section deals with the location and size of the inlet ports. 
Finally, the fifth and last section presents the method used to 
determine the ram weight. 

The pyrotechnic ram device has been designed for normal op- 
eration (separation at 267,000 feet) and its use in an abort sit- 
uation may be somewhat limited due to the excessive drag forces 
on the orbiter which the extended ram members must support. 

These drag forces are the result of the dynamic pressure on the 
orbiter which at nominal separation is only 1 psf. Moreover, 
during the ram operation either both vehicles must be coasting 
or the relative axial velocity of the vehicles must be zero so 
as to minimize the beam bending effect of the ram extension. 

Several constraints were placed upon the physical dimensions 
of the ram and its operation by personnel at MSC. Among the 
most important constraints were the following: 
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1. The output force should vary as the fourth power of the 

elasped time following , pyrotechnic ignitibn, for the 

4 

greater portion of the ram operation. (F = Kt ) 

2. The maximum force should be no greater than 100,000 
lbf . 

3. The operation time should be greater than 1 second, but 
less than 2 seconds. 

4. The stroke should be between 18 and 36 inches. The 
necessary clearance to perform a no recontact separation 
is not yet known. However, the required stroke to 
avoid recontact between the vehicles is probably within 
this range, according to NASA personnel. 



1. Ram Description and Operation 


One of the main advantages of the pyrotechnic ram over the 
other separational techniques investigated is that it could 
most easily be used as the main component of the shuttle mating 
system. Since this investigation is mainly concerned with the 
hardware required to perform vehicle separation, an analysis of 
the suitability of the ram for the mating function will not be 
performed. However, a ram cross section sufficient to withstand 
3.5 times the weight of the orbiter (2.5g), excluding drag, has 
been specified. This is the maximum acceleration which the 
shuttle will experience. 

In Figure (V-l), member 1 provides the ram piston constraint. 
However, it alone will not resist a moment which could rotate one 
vehicle with respect to the other and therefore supports must be 
placed aft to provide this restraint. Moreover, since the shuttle 
is designed for a nominal 2,5 g acceleration and the ram alone, as 
designed, can support only 2.5g, the back supports must take all 
of the longitudinal load caused by the drag on the orbiter. 

Member 2, the ram cylinder, is used to contain member 1 and 
it, along with the bottom section of member 1, serves as the ram 
chamber. The wall thickness of member 2 was determined by con- 
sidering it as a pressure vessel. The maximum pressure it must 
withstand is equal to the output final force divided by the bottom 
cross section of member 1, upon which this pressure must act. 



The magnitude of this pressure for the 6 inch diameter surface 
of member 1 is 3656 psi. 

The function of member 3, the gas delivery chamber, is to 
serve as a manifold to deliver the high pressure gas to the ram 
chamber. Member 3 and member 4 comprise the stagnation chamber. 

This component is welded to member 2. Its wall size is determined 
by considering it as a pressure vessel with the initial pressure 
of 15,000 psi being the maximum that this member must withstand. 

Using this pressure will result in choked flow throughout the ram 
operation. 

Member 4 is the breech where the pyrotechnic is ignited. 

The diameter and length of this vessel was determined by 
considering the density of the pyrotechnic before ignition. 

Member 5, the constraining rod, provides upward transverse 
constraint for the two vehicles throughout the flight. It also stops 
the upward motion of member 1 and thereby concludes the ram operation. 
The frangible nut located on top of this member provides a reliable 
method of allowing the two vehicles to separate. The large spring 
located concentrically around member 5 is used to reduce the large 
tensile shock wave which this component must experience in stopping 
member 1. The smaller spring is used to return member 1 to its 
initial position following separation. The small slot and diameter 
reduction of member 5 are used to provide a means for gas depressuri- 
zation during the ram retraction. 
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Frangible 



PYROTECHNIC RAM, SIDE VIEW, MATED POSITION 


FIGURE V-l 



PYROTECHNIC RAM, FRONT VIEW, FIRED POSITION 


FIGURE V-2 









The material from which the main ram components are fab- 

f 

ricated is Inconel 718. Its former use in a ram on the Apollo 
and its yield strength of 156,000 psi and ultimate tensile 
strength of 183,500 psi at 600°F favor its use in this applica- 
tion. The pyrotechnic used in this analysis is Thiokol TP-Q-3027. 
It has a low flame temperature of 1883°F. The safety factor em- 
ployed in this analysis was 1.5, the same as for the investiga- 
tions of all four separational techniques. The final output 
force for which this was designed is 100,000 pounds. 

In order to commence ram operation, the frangible nut is 
broken by means of two explosive caps contained within the nut. 
Simultaneously, the pyrotechnic material within the breech is 
ignited. The overall time required to reach the maximum pressure 
of 15,000 psi within the stagnation chamber is from 20 to 30 
milleseconds . Initially, gas can enter the ram chamber only 
through the lowest port. As the pressure within this chamber 
builds up, friction is overcome and member 1 proceeds upward. 

As it does, additional inlet ports become available to admit the 
high pressure gas into the chamber. Member 1 is finally stopped 
by the large spring located around member 5. The gas is able to 
exit the ram through the large port at the top of member 1, and 
the smaller return spring returns member 1 to its initial position. 
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2. Separation Dynamics and Calculation of Energy 

Referring to Figure (V-3), the motion of the orbiter is 
described by the following equations. 

Fg 

X ~ c 
o m 

o 

r Fg c 

X * — - dt + c., 

o J m 1 

o 

x o = J( J dt C l) dt + c 2 

o 

where and are integration constants. Since the separation 
force acts through the c.g. of the orbiter, the orbiter motion is 
completely defined by these equations. 

The linear motion of the c.g. of the booster is given by 
Equations (II- 1), (II- 2), and (II- 3). The angular motion of the 
booster can be derived as follows: 


(V-l) 

(V-2) 

(V-3) 



(V-4) 
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Fbg 


dt + c, 


(V-5) 




dt + c 3 ) 


dt + c, 
4 


(V-6) 


where b is as shown in Figure (V-3) 

and and are integration constants. 

The relative velocity between the ram attach points on these 
vehicles is given by 


*■ 4 


X = X + +b 9, 

o DCg b 


(V-7) 


Substituting for X q , , and 0^ from Equations (V-2), 

(II-2), and (V-5), respectively, into Equation (V-7) » leaves 


X * I 


fac J 


F dt + c. 


(V-8) 


where 


i, = r^ + i- + f-i g 

fac L I K J 


and c c is a constant of integration. 
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The relative displacement between the two vehicles is 
given by 


X = X + X, + bQ, 
o beg b 


(V-9) 


Substituting in Equation (V-9) from Equations (V-3), 
(II- 3), an ^ (V-6) for X q , X^, and 0 , respectively, yields 


X = I 


fac 


J(f dt) dt + J 


dt + | c,. dt + Cg 


(V-10) 


where c^ is a constant of integration. In order to determine 
the energy output of this device reference is made to Figure 
(V-4) . The output energy is given by 


E out “ J Fdx + c 7 


(V-ll) 


where c^ is an integration constant. 

For the ramp force from t^ to t^. 


F = K 1 t 


(V-12) 


Output force (Ibf) 
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where is a constant. Subsituting this into Equation (V-10) 
leaves 


X = I 


fac 




t dt ) dt 


) 


(V-13) 


where t^ < t < 

The integration constants c,- and c^ are both equal to zero, 
because the separation velocity and displacement are zero at the 
beginning of motion, respectively. 

Evaluating Equation (V-13) gives 


X = I 


K i t- 


fac 


t 


o 


(V- 14) 


Substituting for F from Equation (V-12) into the above equa- 
tion, yields 




(V-15) 


Therefore, 


F = 


6X 


I £ac t 


(V-X6) 


In order to determine the energy output during this interval, 
Equation (V-ll) is used with F as expressed in Equation (V-16) 
to obtain 


E 


out 2 v 

Ij. t o 
f ac 


r 


X dx 


E 


3 X 


out i, t 2 

fac 


At t = t^ y this becomes 
2 

3X ? 

E 


out 2 

i fac 2 


(V-17) 


(V-18) 


where is the displacement at t = t 

During the interval < t <; t^, the output force is given 

by 

F - K 2 t 4 (V- 19) 

where is a constant. 


The displacement from the beginning of motion may be deter- 


mined by substituting F into Equation (V-10). The integra- 
tion constants, c^ and c^,are the velocity and displacement at 
t = t 2 , respectively. For any time t, in this interval, X is 
obtained as 

X = I fa= C ( f t K 2 t4 dt ) dt 

2 2 



dt + x 2 


(V-20) 


Evaluating and substituting back for F yields 



(t - t 2 )+ X 2 


(V-21) 


from which F becomes 



(V- 22) 


The energy which has gone into the separation at any time 
during this interval may be determined by using Equation (V-ll). 
The integration constant is equal to E 2 , the energy at t = t 2 * 
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Therefore, 


E 


out 



dx + E 2 


30 


I 


IJ- 

facV* 







dx + E 2 



(V-23) 




- + x 2 ] } 


+ E, 


Equation (V-23) gives the total work utilized for separa- 
ting the booster and the orbiter. In order to determine the 
energy which must be used to compare the pyrotechnic ram to the 
other separational techniques we must determine the work done 
against friction. This frictional work is given by 




(V-24) 
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This can be evaluated to yield 


E 


fr 


F ~ X Q 
fr 3 


(V-25) 


The energy capability of this device is given by 


E = 



+ E 


out 


(V-26) 


In utilizing the preceding equations for the energy output 
and the vehicles’ relative dynamics, the variables and K 2 
are defined as 



(V-27) 



(V-28) 


Another variable of use is the time factor, which specifies the 
ratio of the total amount of time during which the ram operates 
be the time during which the ramp portion is in operation. 


fc fac t„ 


(V-29) 
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Finally, in order that the ramp function force and the fourth 
degree function force be equal at t^, the following relation 
must hold: 


K 1 t 2 " K 2 l 2 


(V-30) 


Among these four equations, (V-27), (V-28), (V-29), and (V-30), 
there are seven variables. If the quantities t£ ac > and 
are specified, they may be solved simultaneously. 

For the thermodynamics involved in the ram operation, the 
work involved in moving the ram pisLon from one point to the 
next must be determined. The differential work going to fric- 
tion must also be calculated for each point. 

Therefore, for the interval t^ < t c!^, the equation is 


A Work. 


JC. *x. r x. 

. = I 3 F. dx - j J F. - dx + I J F £ 

J J i J j.-l fr 

w o ■ o J X., 

t-1 


dx 


(V-31) 


where 


•j .ref erg to any iteration step, 

j“l refers to the preceding iteration' step, 


F 


fr 


is the frictional force, 


F & F - must both be evaluated using Equation (V-16) 

j j' 1 


Substituting into Equation (V-31) for F . and F from 

J J 

Equation (V-16), and integrating, leaves 


2 2 
X . X 


&WOrk j = Ffr <X j ' X M } 


(V-32) 


j J“1 


During the interval t^ < t < t^, the differential work is 
given by 


r x. X X. 

AWork. = J F.dx- J F . - dx + | J F dx (V-33) 

J J J y 

Art Aa A , i 


x m 


Evaluating F and by using Equation (V-22) and sub- 

stituting, leaves 


AWork. = 


30 


j I 


f ac 


{ (t 2 1 . f jQ [ X - X 2 (tj - V - X 2 ] dx + 


(t 


X 

[ x - x 2 (‘j.1 - V - x 2 ] dx } 

j-1 “ fc 2 > X 2 


I 


+ 3 F^ dx 


x j-i 


fr 


(V-34) 


Integrating and evaluating Equation (V-34) leaves 


A Work . 


J 



- < X j - X 2> (X 2 (t j ' t 2> + V ] + 



-CX..1- x 2 ) (x 2 (t._ 1 - t 2 ) + x 2 


>]} 


+ F fr (X j - X j-1> 


(V-35) 
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3. Operational Thermodynamics 

The method used to size and locate the inlet ports requires 
the calculation of the pressure, temperature, and mass of the 

1 

: f‘ pyrotechnic medium in the ram chamber as shown in Figure (V-5). 

Several assumptions had to be made in order to analyze the 
thermodynamics. First, it was necessary to assume that the pryo- 
technic medium is an ideal gas. This assumption is approximately 
true for the majority of the constituent gases. If probably does 
not introduce too great of an error into the calculations. 

This assumption implies that the Equation of State 


Pv = RT (V-36) 

can be applied. It also means that the internal energy and the 
enthalpy are functions of temperature only so that they are governed 
by the equations 


u = C T 
v 


and 


(V-37) 


h = C T 
P 
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FIGURE V-5 




Second, the process is essentially adiabatic. Since the 
only means of heat transfer from the control volume is by conduc- 
tion through the ram walls and since the process occurs during a 
very short (1.0 - 2.0 sec) time interval this assumption is valid. 

Third, the losses due to friction are essentially zero. 

This assumption is reasonable because the losses due to frictional 
flow are very small except for the first inlet port and the firc- 
tional energy loss due to the movement of the ram piston is only 
1.07o of the total. This assumption, along with assumption two, 
imply that the process is isentropic so that the equation 

Pv k = C (V-38) 

where C is a constant, applies. 

In addition to these assumptions the constraint is put on 
the process that only choked flow occurs across the inlet ports. 

According to Reference (3), the First Law of Thermodynamics 
for an unsteady process occurring during a small time interval, 

A t, is given by 



Q + S m 

c.v. 


i ( h i 


Vel 


g. 


+ Z 


&J 
i 8, 


)- 


Work + 

c.v. 


+ s m e 



+ 



+ Z e f-) + [»! («*! 

c 


■f 



+ 



- m ( u 
o \ 


Vel 


o + 2 g. 


+ Z 


o g. 


:)] 


c.v, 


(V-39) 


where 

the subscript i refers to quantities entering the control 

volume during ^t 

e refers to quantities leaving the control 
volume during 

1 refers to quantities within the control 
volume after £t 

0 refers to quantities within the control 
volume before At 

Q is the heat flowing out of the control 

C • v • 

volume, 

Work is the work output of the process 

. C • v • 

Z's refer to the potential energy associated 
with changes of elevation. 
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In applying Equation (V-39) across the ram inlet ports, 
Figure (V-5), the following conditions are considered. 

(a) Q =0 because of the adiabatic assumption, 

C • V • 

(b) The change in elevation is negligible, Z ^=Z ^=Z ^=Z ^=0 . r 

(c) The mass exiting from the expanding ram chamber is 
zero, 2m =0. (until the ram motion is complete), 

(d) The velocity of the gas within the ram chamber is 

negligible so that Vel. = Vel = 0. 

1 o 

With these considerations. Equation (V-39) reduces to 

2 

. Vel. v 

2 m. (h. 4* — — — ) = Work + fm. u. - m u 1 (V-40) 

i \ i 2 g / c.v. L 1 1 ooJ 


Substituting for h^, and u^ from Equation (V-38) gives 
/ Vel A - 

2 m. (C T. + - — — ) = Work + m.C T - m C T (V-41) 
i\pi 2 g / c.v. L 1 v 1 ovoJ 

c 

Equation (V-41) describes the general problem of gas entering an 
adiabatic system with movable boundaries, but with no gas escaping 
from this expandable control volume. Using Figure (V-5), Equation 
(V-41) may be rewritten as 

2 


2 m 


t, 1 (' 


Vel. 


C T + -~-=- 
p t,0 2 g 


^ = AWork^ + 


(V-42) 


[ m r,l 


CT . - m _ C T 
v r,l r,0 v r 


j 


where 


Z m 


t,l 


t,o 


Vel 


t,0 


^Work. 


m 


r,l 


r,l 


m 


r,0 


r,0 


= summation of the mass of gas which has entered 
in the port throats at step number 1, lbm 
= temperature of the gas in the port throats at 
step number 0, °R 

= velocity of the gas in the port throats at 
step number 0, in/s 

= differential work in moving the ram piston, 
lbf-in 


= the mass of gas in the ram chamber at step 
number 1, lbm 

= the temperature of the gas in the ram chamber 
at step number 1, °R 

= the mass of gas in the ram chamber at step 
number 0, lbm 

= the temperature of the gas in the ram chamber 
at step number 0, °R 

Step 0 corresponds to the time when no gas has yet entered 
the ram chamber. Step! corresponds to the time when the fric- 
tion opposing the. motion of the ram piston has just been over- 


come. 


Since at step 0, no gas is in the ram chamber, m^ ^ = 0. 
Because the ram piston has not moved 9 friction has not 
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been overcome until right at Step 1, and therefore AWork^ = 0. 
Equation (V-42) therefore reduces to 


S m 


t,l (' 


Vel 


C T A + 
P t,0 


2 8 , 


t,0 




m - C T - 
r, 1 v r, 1 


(V-43) 


From the continuity equation, however. 




(V-44) 


This means that the summation of the mass of gas which has 
entered through the ports at Step 1 is equal to the mass of gas 
in the ram chamber at Step 1. Furthermore, since the tempera- 
ture, q, and the velocity, Vel fc q, in the throat of each of 
these ports is the same, we can obtain the following expression 
for the ram chamber gas temperature immediately preceding the 
onset of motion. 



Vel t 0 

l t,0 + 2C g 


VC 


(V-45) 


where ]£ is the specific heat ratio, C /C 

p v 

Since both T t q and Vel t Q are decreasing with respect to 
time. Equation (V-45) is valid only for a small time interval 
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At. In order to determine T. A and Vel. n > Reference (4) gives 

t , 0 t , 0 


T t,0 ■ < T fac>( T s,0> 
Vel t,0 - V s c R T t,0 


(V-46) 


where T^^is the ratio of choked temperature to stagnation tem- 
perature, commonly called the temperature factor. 

After some gas has entered the ram chamber. Equation (V-42) 
may then be used to determine the gas properties. Writing it 
in a general form yields 



+ 



T. . , + 

t,j -l 


C 

v 


T , 
r O 


Vel t j-i \ 

-Ft 7-) = 

“r.j -1- °v T r, j-ll 


(V-47) 


where j. is the step number. 

In this equation there will, in general, be three unknowns, 
2 m , m and T However, from the continuity equation 

one may conclude that 


E m 




= m 


- m 


.j r « 3-1 


(V-48) 


Applying the Equation of State leaves 




P ,Vol , 
r ».1 * V 1 

R 


(V-49) 


where P . and Vol . can be evaluated by knowing the force on 
r » J r > J 

the ram at step j and the relative displacement at j, respec- 
tively. 

Substituting for Em from Equation (V-48) and for m 

*- > 1 r s j 

from Equation (V-49) into Equation (V- 47), one can obtain the 
mass within the ram chamber as 


2 S, 


m . = 

r,j 


2g C T . - + Vel . . L J 
c p t,j-l t,j-l 


P Vol . 

0 -{AWork. + C 4- 

2 l i R v 


Vel^ , 


+ m (CT - CT - + — 

r,j-l \ p t,j-l v r,j-l 2g )] 


(V-50) 


The gas temperature corresponding to this mass is 


P , Vol . 

T a- .1 J 

r >j m .R 
r, j 


(V-51) 


The mass of gas remaining in the stagnation chamber is given by 


the continuity equation. 
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m , = m - m . 
s,j s,0 r,j 


(V-52) 


where s refers to the stagnation changer and m n is the mass of 

S 9 V/ 

gas within the stagnation chamber at step 0. It may be determined 
by substituting the pressure, temperature, and specific volume of 
this gas at step 0 into the Equation at State. 

The general equation for the pressure in the stagnation chamber is 


,Vol ,-k 


P . = C 

s,J 


e) 

S,J 


(V-53) 


where C is evaluated before any gas has entered the ram chamber. 

The temperature in the stagnation chamber can be calculated 
from the Equation of State: 


P . Vol 
■id S 


s,J 


R m 


(V-54) 




Finally, the throat temperature in the ports is given by 


T t,j = 


(V-55) 


and the gas velocity, which is at Mach 1, is given by 


Vel . « k g RT . 

t,j 1 & c t ,3 


(V-56) 
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These two quantities T . and Vel . can be used in Equa- 

8 s j J 

tion (V-50) to determine the gas properties in the next iteration 

step* 

The greatest source of error in the calculations for the 
ram operation centers around the assumptions made at the be- 
ginning of this section. Additionally, the small diameters of 
the first few ports causes the frictional losses accompanying 
flow through a narrow port to become dominant. For that reason 
the length to throat diameter ratio for the smallest size inlet 
port (port 1) was not allowed to be greater than 20. This pro- 
duces about a 157o deviation between the calculated and actual 
pressure drop across this port and about a 12% deviation between 
the calculated and actual mass flow rate. By fabricating the 
ports in the form of converging-diverging nozzles, this error 
may be reduced since the length to throat diameter ratio is 
smaller. It is recommended that some experimental work be per- 
formed on this design in order to compensate for the slight dis- 
crepancies which exist between the actual gas and its flow 
characteristics and the idealized assumptions which had to be 
made in order to mathematical lysiimulate the ram operation. 


0 - 
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4. Port Sizes and Locations 

The size of the first port is assumed because it is a free 
variable. If it is assumed too large the operation time will be 
less than 1 second. If it is assumed too small the operation 
time becomes greater than 2 seconds. These conditions are 
contrary to the time constraint explained at the beginning of 
this chapter. The location of the first port is as shown in 
Figure (V-5). The initial mass flow rate into the ram chamber 
is only through this point. According to Reference (5) it is 
given by 



s, 1 


— A s 

1 l FT 

i s, 1 


m 


fac * 


(V-57) 


where A. is the area of port 1 and 




1 

k+1 (k + ])/2(k - 1) 

<-y) 


Since m , can be easily evaluated from the Equation of State, 
(V-49) , the time required to build up enough pressure in the ram 
chamber in order to overcome friction is given by 





(V-58) 
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For any interval after interval 1, the general equation for 
the mass flowing into the ram chamber is given by 



(V-59) 


(V- 60) 


the increase in total port area which must occur in order that 
the force output function be maintained can be calculated. 

The mass flow rate of any step in the ram operation may be 
evaluated as 


m 


r,j 


m . , - m . , 
r a ,i+l r ,.1-1 

2At 


(V-61) 


Substitution of Equation (V-60) into Equation (V-59) and 
rearranging gives 



m . \ T , 

Lxll 

V* m fac 



(V-62) 



The location of any port is at the displacement corresponding 

to m . , T . , and P . At certain points during the ram 
s,j s,j 

operation, evaluation of Equation (V-62) will yield a negative 
value for . This is because this position does not require 
the location of an inlet port. That is, more high pressure gas 
is not required to maintain the desired output force function. 

As can be seen from Appendix B, these points are only encountered 
during the time from t^ until t^ . It is felt that since the 
amount of energy which the ram puts into the system during this 
ramp portion is only a negligible portion of the total, the 
deviation from actual operation is not significant. Again, as 
explained at the end of the thermodynamic section of this chapter, 
the computed port sizes represent the lower limit of the actual 
sizes required . This is due to the assumptions of isentropic 
flow of an ideal gas. 


5. Ram Weight 


The method used in determining ram weight consisted in 
weighing the ram shown in Figure (V-6). This ram contains the 
predominant components of the ram shown in Figures (V-l) and 
(V- 2) . The terminology used in the computer program shown in 
Appendix A and Appendix B to determine the ram size and weight 
is illustrated in Figure (V-6). The resulting dimensions of a 
ram capable of supplying 130,207 ft-lbf to the booster and orbiter 
separation is shown in Figure (V-7). In this approach the weight 
of the stop spring, return spring, explosive nut, required washers, 
piston rings, seals, and weld material are compensated for by 
multiplying the resultant weight by an allowance factor of 1.20. 

It is felt that this 20% addition to the ram shown in Figure (V-6) 
is a conservative approximation because certain aspects of the 
model made it heavier, such as the reduction in member 5 area 
and the slight effect of the ports. 

The ram weight was determined by knowing the ram stroke, 
either an inside or outside diameter for each member, the design 
pressures, the safety factor, and the metal tensile strength and 
density. The ram stroke determined the lengths of each member and 
the thicknesses of the pressure vessels could be determined from 
knowing the inside radius, pressure, and the tensile strength. 

The volume of each member was computed and multiplied by the 
density to obtain the weight. 









CHAPTER VI 


RESULTS 

The purpose of this study was to determine the minimum weight 
separation technique among the four investigated. Since the 
energy required for separation is a function of the final 
translational velocities of both vehicles and since these will not 
be known until the shuttle design is further refined, it is much 
more useful to show the graphical relationship between weight and 
energy than it would be to assume final specific translation 
velocities. Figures (VI- 1) and (VI-2) are the results of this 
study. Figure (VI- 1) shows all four separation systems and 
Figure (VI- 2) decreases the range of the abscissa in order to 
better display the two least weight techniques. A safety factor 
of 1.5 is used for each technique because the safety factor 
employed in the design of the rocket engine cases used in the 
analysis of the rocket separational scheme was 1.5. The excessive 
weights of the spring for the spring loaded ram and the ejector 
gas in the case of the high speed, non-ignitable gas ejectors 
rules out their further study. The data for the rocket engines 
follows a range of values depending upon which numbers are taken 
as the vacuum specific impulse and the mass fraction. The curve 
for the pyrotechnic ram reveals that it is only good for a range 
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of values from 80,000 ft-lbf to 160,000 ft-lbf. This is due to 
the four constraints placed on the ram and its operation by MSC 
personnel, explained at the beginning of Chapter V. 

Referring to Figure (VI- 2) it may be noted that the use of 
rockets results in a lighter separation subsystem design but for 
the range of values for which the ram may be used, it would result 
in a lighter overall shuttle system design since it would be dual 
functional. Because the method used to constrain the vehicles 
together before separation will employ one forward constraint and 
a set of aft supports, the ram piston could be conveniently used 
for the forward constraint. The pyrotechnic technique is the 
only separation method investigated which has this capability. If 
it can be said that 45% or greater of the ram weight is the weight 
which would be required to mate these two vehicles, then it can be 
shown that the pyrotechnic ram would be the most feasible 
separation technique. 
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APPENDIX A 


COMPUTER PROGRAM LISTING" 





** 




COMPUTER FLOWCHART 







COMMENT ********** PYROTECHNIC RAM PROGRAM ********** 

COMMENT ********** DATA INPUT ********** 

( 

IMPLICIT REAL <L,M,K) 

REAL IB,INRFAC,INRSYS 

DIMENSION AREACI25) , AREAT ( 25 ) , DI AT ( 25 ) , FOUT ( 25 ) » FR AM< 25 ) , 
1FRAMPI25) ,F4DEG<25) ,MASSDR(25) , MASSRt 25 ) , MASSS ( 25 ) , PRESR ( 25 ) , 
2PRESSC25) , RHOS (251 , RHQT ( 25 ) , TE MPR( 25)/TEMPS(25)» TEMPT 1251* 
3TIMEI 25) , TIME AM ( 25 ) , TI MEAA( 25 ) , T I MEFFI 251 , VOLRC 25) f 
4VELGT (25) , WORKDL (25),X(25),XD( 25),NUM< 251 
READ! 5 1 1) KGtMOLt TEMFLMf SPHTPR, SPHTVL 
READ! 5,1 } CRAMP, FFR,FOUTF,T I MPER 
READ( 5,1) THDF AC , VOLFAC 
READ! 5*2) B, IB, MB, MO 
READ! 5,3) NSTEPS 

READ! 5,1) PRES SI # RHO,SFTFAC, STRSSU, STRSSY, WTEF AC 
READ( 5,1)SPACE1 

READ(5,1)DIARM,LN11,LTSFAC,LT1FAC,THC11 

READ(5,1)D23FAC,LT2FAC 

READ! 5 , 1 ) LN32 , LT3F AC 

READ! 5,1) DIA^l ,LN41*LTAFAC 

READ(5,1)DIA52,LN51,LN52 

1 FORMAT ( 8F10« 3 ) 

2 FORMAT ( F10.3,7F10*C) 

3 FORMAT (110) 


ON 

Ln 


COMMENT ********** EVALUATING CONSTANTS ********** 


P I =3c 14159 

NSTEPM=NSTEPS-1 

ASTEPS=NSTEPS 

K1=1/KG 

K 2 = ( K G - 1 * }*K1 

K3=l • /K2 

K 4 = 1 . /« l.-KG) 

K5=-K4 

TEMPSI=TEMFLM+460* 

R=1545.*12®/M0L 

G2J 144=32. 17*2*778. *144. 

CRTKG1-SQRT(KG*32. 17*12/R) 

CRTKG2=(KG+1. »/2. 

CRTKG 3= < KG+1. ) /<2*(KG-1. H 
MDFAC =CRTKSl*(l./{ {CRTKG2 ) **CRTKG3 ) ) 
RHOFAC= ( 2/ ( KG+ 1. ) )**K5 
PRSFAC=(2/(KG+ls ) )**K3 

tmpfac-?2/(kg+i. n 

GC=32« 17 
IB=IB*GC 

INRFAC=GC*12.*U 1. /MO) + U . /MB) +< B**2/IB» ) 

INRSYS=12./INRFAC 

TIMPER=TIMPER/1G0. 

timfac=i./timper 

frampf=foutf/timfac**4 

timea=frampf/cramp 

timef=timea/timper 

timper=timper*igo. 

timfac=tinfac-.oi 

A J A = ( (ASTEPS-1. )/TIMFAC)+l. 

timfac=timfac+.gi 

J A=A J A 

C4DEG=F0UTF/TIMEF**4 



XA=INRFAC*CRAMP*TIMEA**3/6. 
X0A=INRFAC*CRAMP*TIMEA**2/2, 
ENERGA=3 0 * : XA* :5 * e 2/(INRFAC 3 ! c TIMEA^*2*12* ) 

XF = ( I NRFAC *C4D-EG* (TIMEF**6-TIMEA*’*'6)/30* ) 
1+XDA*(TIMEF-TIMEA)+XA 


COMMENT ********* RAM 


SIZING AND WEIGHING ********** 


AREARM=(PI/4. I *{ DI ARM**2-D I A52 **2 ) 

DI Al2=DI ARM 

V0LRI=SPACEl*AREARM/4. 
VOLRF=AREARM*XF+VOLRI 
VOLS= VOLRF*VOLFAC 
STR SSD=STR SSU/ SFTF AC 
LN12=LT$FAC*XF 
LN31 =XF 

LN30=LN31+LN32 

VOLS2=( PI*DIA41**2/4. i*LN41 

VOLSl=VOLS-VOLS2 

DIA31=SQRT(4.*VOLS1/(PI*LN30 H 

PRESRF=(FOUTF+FFR)/AREARM 

DI A2 1=DI A12 

LN21=LN12+SPACE1 

THC21=(PRESRF*DIA21>/(2.*STRSSD> 

THC22=THC21*LT2FAC 

DIA22=DIA2U-2.*THC21 

DIA23=DIA22*D23FAC 

THC31=(PRE$SI*DIA31>/(2.*S7RSSD) 

THC32=THC31*LT3FAC 

THC41=(PRESSI*DIA41)/(2.*STRSSD| 

THC42-THC4I*LT4FAC 

DI A42 -DIA41+2* *THC41 

DIA32=DIA31+2.*THC31 

THC12=THCI1*LT1FAC 

LN1 0=LN1 1+LN1 2 

LN50 = LN21+LN5H-LN11+LN52 

DIA1I=DIAI2“2.*THCI1 

V0L11=(PI*(DIA12**2-DIA11**2 )/4, )*LNlO 
VQL1 2 = ( PI *(DI A 11**2*DI A52**2 )/4. )*THC12*2. 
VOL 1= VOL l 1 + VOL 12 

V0L2i-(PI*(DIA22**2-DIA21**2 >/4. )*LN21 


VOL22=(PI*(DI A23**2-DI A52**2 >/4. )*THC22 
V0L2=V0L21+V0L22 

V0L31=(PI*(DI A32**2-DI A3 1**2 1/4, )*LN3C 
V0L3=V0L31 

V0L41=(PI*(DI A42**2-DIA41**2 )/4. )*LN41 

V0L42=(PI*(DIA42**2 1/4. )*THC42*2. 

V0L4=V0L41+V0L42 

V0L5=(PI*0IA52**2/4. )*LN50 

WT1*RHO*V0L1 

WT2=RHO*VOL2 

WT3-RH0%V0L3 

WT4-RH0+V0L4 

WT5=RHG*V0L5 

WTMEM=WT1+WT2+WT3+WT4+WT5 
WTERR=WTMEM*WTEFAC 
MASSSI=PRESSl*VOLS /<TEMPSI*R> 

WTGAS=MASSSI 

WEIGHT=WTMEM+WTGAS*WTERR 
0 1 AMI N=THC2 1/ THOFAC 
ARAM I N= PI I AM IN**2/4. 

ENERGF= ( ( 30. / ( INRFAC* I T I MEF**2-T IME A**2 I ) ) 
l*t ( <XF**2-XA**2)/2. )-( XF-XA) *( XOA* 

2 ( TI MEF-TIMEA ) + XA ) ) J/12.+ENERGA 



COMMENT ********** GAS CONDITIONS 


AT D AND 1 ********** 


CQNGAS=PRESSI*<MASSSI/VOLS)**(-KG) 

PRESSF-PRESRF/ PRSF AC 

TEMPTI=TMPFAC*TEMPS I 

VELGTI=SQRT(KG*32. 17*12. *R*T EMPT I) 

TEMPR(1)=KG*TEMPTI*(VELGTI**2/(G2J1A4*SPHTVU ) 

VOLR ( 1 ) =VOLRI 

PRE SR ( 1 ) = FFR/ ARE ARM 

MASSRC 1 ) =PRESR( 1 1 *VOLR (1 ) / ( R*TEMPR( 1 ) ) 

MASSSC 1 )=MASSS I-MASSRI 1) 

PRESS < 1 ) =CONG AS* (M ASSS { 1 ) /VOLS ) **KG 
TEMPS (1 IMPRESS Cl)*VOLS/(R*MASSS( 1 ) ) 

RHOSm = < PRESS m /CONGAS >**K1 
RHOT ( 1 ) -RHOS ( 1 ) *RHOF AC 
AREAT(1)=ARA,MIN 

DT AT( 1 ) -SORT ( ABS < AREAT( 1 »*4. /PI)) 

TEMPT( 1 J^TEMPS ( 1 ) *TMPFAC 

VELGT(I) = SQRT(KG*GC*12.*R*TEMPT( 1) ) 

PSLSQT -PRE SS ( 1 1 /SQRTC TEMPS ( 1 ) ) 

MAS SDR (1)=AREAT(1 )*( PSLSQT *MDF AC ) 

TIMEBU=MASSR( l)/MASSDR(l) 

TlMOPR=TIMEF+T IMEBU 


COMMENT ***** DYNAMICStPRESSURE? VOLJME, AND DELTA WORK GENERATION ***** 


COMMENT ********** FORCE OUTPUT IS A RAMP FUNCTION ********** 

TMINCR=TIMEF/( ASTEPS-1. ) 

X<11=0® 

XD(n=0. 

FRAMf I } =FFR 
TIMEAAC 1)=0« 

DO 400J=2,JA 

TIMEAAI J) = TIMEAA( J-D+TMINCR 
TIMEAMI J)=TI MEAAC J 1 
T I ME( J)~TIME AMI J) 4-TIME BU 
FRAMPU)=CRAMP*TIMEAA( J) 

FOUT ( J )=FRAMPC J ) 

FRAMI J)=FOUT« JI+FFR 

XDC J ) = 1 NRFAC*C RAMP*TI MEAAC J I **2/2» 

XC J)=INRFAC*CRAMP*TIMEAA( Ji**3/6. 

PRE SR ( j ) = FRAM( JS/AREARM 
VOLRC J* =C X< J| *AREARM) +VOLRI 
IFCJ#NE®2) GO TO 402 

WORKDLI J) = <3./INRFAC)*C (XC J ) **2/TI MEAA ( J >**2l) 
i+FFR*X{ J) 

GO TO 400 

402 W0RKDL(J) = (3*/INRFAC)*((X(J$**2/TIMEAA( J 1**21 
1-<X< J~1)**2/TIMEAA( J-1)**2H+FFR*(XI J )-X< J-l ) ) 

400 CONTINUE 


COMMENT ********** FORCE OUTPUT IS a 4 DEGREE FUNCTION ********** 


J A1 = J A+ 1 

TIMEFF ( JA)=TIMEAA( J Al 
DO 401 J=JA1,NSTEPS 
A J = J- 1 

TIMEFF ( J) =Tl MEP*A JV ( ASTEPS — 1 * ) 

T I ME A M ( J ) = T I M E F F ( J ) 

TIME { J)=TlMEAM< J)+TIMEBU 
F4DEG(J)=C4DEG*TIMEFF(J) **4 

XD ( J ) = ( INRFAC*C4DEG*( T I MEFF( J) **5-T I MEA**5 ) /5m ) 

1 + XD ( J A) 

X< J) = (I NRFAC*C4DEG*(TIMEFF< J)**6-TIMEA**6l/30. ) 
l + XDI JA)*<TIMEFF< J)-TIMEA)*X< JA ) 

FOUT( J )=F4DEG ( J ) 

FRAM1 J ) =FOUT ( J ) +FFR 
PRE SR ( J )=FRAM ( JJ/AREARM 
VOLR < J) =X ( J) *AREARM+VOLRi' 

WORKDLl J> =f 3D. /INRFAC)*( (l./CTIMEFFf J)**2-TIMEA 
1**2) ) *( <(X( J)**2-X( JA)**2)/2»)-tXf J)-X< JA) )* 
2(XD< J A) *< TIMEFF <J ) -Tl MEA ) +X ( JA ) ) )-< l./CTIMEFFU- 
31)**2-TIMEA**2) )*( ( (X ( J-l I<**2-X( JA) **2 )/2* ) 

4- ( X( J-l )-X < J A I ) * ( XD ( JA) * ( T I MEFF( J-l ) -TIME A) 

5+X( JA) )))+FFR*<X( J)-X(J-l) ) 

401 CONTINUE 

TIME ( 1 1 — T I MEBU 
TI MEAM( 1 ) =0* 

FOUT ( 1 ) =0» 



COMMENT ********** THFRMODYNAMI CS GENERATION ********** 
DO 500 J= 2 1 NSTEPS 

MSDUM1=S-PHTVL*UPRESR< J**VOLR( J) /R> -MASSR( J-ll* 

IT EM PR ( J-l I ) 

MSDUM2=MASSR( J-l )* ( SPHTPR*TE MPT ( J-l ) + ( VELGT( J-l ) 
1**2/G2J144> ) 

MSDUM3-MSDUM2/MASSR( J-l ) 

MASSR(J)=( ( WORKDL (J)/(12«*778* ) ) +MSDUM1+MSD JM2 ) / 
1MSDUM3 

TEMPR ( J ) = PRESR ( J )*VOLR ( J > / C MASSR( J ) *R) 

MASSS(J)- MASS'S I- MASSR(J) 

PRESS ( J) = CONGAS*<MASS$U) /VOLS )**KG 
TEMPS! J IMPRESS < J ) *VOLS./ <R*MASSS< J I ) 

RHOSi J) =( PRESS ( J ) /CONGAS I **K1 

RHOT C J ) =RHOS ( J ) +RHOFAC 

TEMPT ( J )- TEMPS ( J) *TMPFAC 

VELGT ( JMSQRT! KG*32* 17*12* *R*T EM PT( J ) ) 

500 CONTINUE 


COMMENT ********** PORT SIZES AND LOCATIONS ********** 

AMI SUM= ARE AT ( I ) 

NUM ( 1 ) = 1 
NSUM= 1 

AREAC ( 1 ) = ARFAT ( 1 ) 

DO 501 J=2 » NSTEPS 
PSL$QT=PRESS( j) /SORT (TEMPS ( J )) 

I F( J-NSTE PM) 503, 503,504 

503 MASS DR ( J ) = ( MASSR ( J+1)-MA$SR( J-l) ) / ( 2 . *TM INCR ) 

GO TO 505 

504 MAS SDR ( J)=(MASSR(J ) -MAS SR (J-l) )/TMINCR 

50 5 AREAT ( J ) = ( MASSDRC J ) / ( PSLSQT*MDFAC ) )-AMlSUM 
IF ( AREAK J ) ) 506 9 506,507 

506 AREAT(J)~ 0*000 000 
DIATC J) =0*003000 
ARE AC (J)=AREAC(J-1) 

NUM( J ) =0 

GO TO 501 

507 IF( AREAK J) -ARAM IN) 506, 506,508 

508 AMI SUM-AMI SUM* ARE AT ( J ) 

DlATC J)=SQRT(4 # *AREAT( J)/PI ) 

NSUM=NSUM+ 1 

NUM( J)=NSUM 
AREAC( J)=AMlSUM 
501 CONTINUE 


* 


l 


t 


COMMENT **********P R I NTOUT ********** PRINTOUT ********** 
COMMENT ********** PRINTOUT IMPORTANT PARAMETERS ********** 


WRITE ( 6, 1C ) 

ID FORMAT! 1H1 ,5! /) ,T52,' PYROTECHNIC RAM ANALYSIS*//// /) 

WRITE (6tll)ENERGF, WEIGHT 

II FORMAT! T48, 'DATA POINT-ENERGY VS WEIGHT CURVE* // 

1T49 , * ENERGY (FT-LBF ) * * 

2T70,*WEIGHT(LB)*/ T50 ,F10. 1 , T67, FlO« I ////) 

WRITE (6»2G )B»IBr MB » MO,INRSYS 

20 FORMAT ( T53, * CONFIGURATION P AR AMETERS • ///T4 , • DIST.B! FT) # * 

l T 24 , • I NER* BOOSTER! LBM- FT* *2)* ,T56 f • MASS BOOSTER! L BM )• , 

2T83 , * MASS ORB ITER! LBM) • ,T108 INER. SYS.(SLUG)'/ 

3T2 ,Fl0.3,T26,F15.0 f T53 f Fl5.0,T8G,F15.0,Tl04,F15.0////) 

WR I TE ( 6 ,30 ) 

30 FORMAT! T5C » * R AM OPERATIONAL CHARACTERISTICS 1 // > 
WRITE(6,31)TIMPER,TIM0PR,TMINCR, T I ME BU ,Tl MEA, T I MEF 

31 FORMAT! T7 , • TI ME PERC. RAMP' , T27, • T I ME OPER.(SEC)* 

1 T47* 'TIME INCR. (SEC)' ,T68, 'TIME l ( SEC ) ' , T89 , • T I ME 2(SEC)', 
2T109, 'TIME 3<SEC>'/ 

3T7,F10.2,T27,F10.4,T47,F10.4,T67,F10.4,T88,F10.4,T1G8,F10.4// ) 
ENERG2-ENERGA + !FFR*XA/12# ) 
WRITE(6,32)CRAMP,FRAMPF,XA,XDA,ENERG2 

32 FORMAT! T8,*K1 (LBF/SEC) • f T33, 9 FORCE 2( LBF ) ' , 

1T60,'DISP.2< IN) ' ,T84, 'VEL.2! IN/SEC) * , T 10 7, * ENERGY 2! LBF-FT ) * / 
2T9,F1G.2,T33, F ID, 2 , T59 , FI 0. 5 ,T 84, F10. 5,T106, F12. 2// ) 

ENERG3=ENERGF*(FER*XF/12. ) 

WRITE (6,33) C4D EG, FOUTF»XFt XD(N ST EPS ) , ENERG3 

33 FORMAT! T6,*K2 ( L8F/SEC**4 )* ,T 33 ,* FORCE 3! LBF ) ' , 

1T60, 'STROKE! IN)* , T84, • VEL. 3 ( IN/S EC ) ' , T 107 ,* ENERGY 3! LBF-FT) •/ 
2T7,F10* 2,T33,F10.2,T59,F10. 5,T85, F10.5,T108,F10.2//) 
WRITE(6,34)FFR,N STEPS 

34 FORMAT! T40,'FRIC. FORCE ( LBF )' ,T68, 'NUMBER UF STEPS ANALYZED'/ 
1 T40 * F 1 C • 2 , T76,I5////) 


--j 

Ln 



n 


WRITE<6,40> 

40 FORMAT ( 1H 1 , 6 { / ) , T46 , • PYROTECH NI C GAS ( TP-Q-30 27 ) PROP ER I T I ES' // ) 
WRITE(6,41 >MOL,TEMFLM,R,KG,SPHTPR,SPHTVL 

41 FORMAT ( T5 , » MOL« WE IGHT • ,T22 t * FLAME TEMP. ( F ) » , T4 1 t 

1 * GAS CONST. -R (IN-LBF/LBM-R) ' ,T72, 'SP. HEAT RATIO' ,T9l, 

2 'CPOC BTU/L8M-R ) * ,T110 ,' CVO ( BTJ/LBM-R ) 9 / 

3T4,F10.5,T21,F10.C,T47,F19.2 ,T 74, F8. 4 ,T93 , F8. 4, T 1 1 2, F8. 4/ / I 
WRITE (6 ,42 )RHOFAC, PRSF AC , TMPFAC , CONGAS » MDF AC 

42 FORMAT ( T6 ,'OEN. FACTOR' , T30, * PRES. FACTOR', T54, 

1 • TEMP. FACTOR' ,T80, ’STAG. PRES. FACTOR' ,T112, • FLOW FACTOR' / 

2T5,F1G.4,T29,F10.4,T53,F10.4,T81,F12.0,T110,F10.4 //// ) 

WRI TE ( 6, 5G ) 

50 FORMAT (- T55,'RAM DESIGN PARAMETERS' // ) 

WRITE(6,5nPRESRF,PRESSI ,SFTFAC 

51 FORMAT ( T1 6, 'DESIGN PRES.-MEM.2(LBF/IN**2)* » T6 1, 

1 'DESIGN PRES. -HEMS. 3 AND 4( LBF/ IN**2 )' , T58, 'SAFETY FACTOR'/ 

2T23,F1G.1,T92, F10.lt T56tF10.2//) 

WRITE(6,52)STRSSU,STRSSY,RH0 

52 FORMAT ( TIOt'RAM MATERIAL' ,Y32, • ULTIMATE STRENGTH! LBF/IN**2 ) • ,T70, 

1' YIELD P0INT(LBF/IN**2I' ,T103, ' OENS ITY ( LBM/ I N**3) ' / 

2T10, 'INCONEL 718' , T40 , F10. 0 , T74, F10. 0 , T 104, FI 0. 4// // ) 

WR I TE C 6,60) 

60 FORMAT! T49,'RAM SIZING AND OPERATIONAL FACTORS'// ) 

WRITE (6,61 ) VOL FAC ,WTEF AC, THDFAC, INRF ACtTI MFAC 

61 FORMAT! T9 ,' VOLF AC' ,T29 ,' ERROR FACTOR' , T54, ' THC/PQRT DIA.' 

1T79, • INERTIA F ACTOR' , T107 , » TIMEF/TI ME A* / 

2T9,F5.2 ,T32,F5.2,T58,F5.2,T80,F10.6,T1G5*F9.2 // ) 

WRITE (6, 62) LT1FAC,LT2FAC,LT3FAC,LT4FAC,LT$FAC,D23FAC 

62 FORMAT ( Til , ’ LT1FAC' ,T31 , ' LT2FAC* ,T51» 'LT3FAC* ,T71, ' LT4FAC ' , 

1 T91 , ' LTSFAC « , T 1 1 1 , * D23FAC • / 

2T7,F9.2,T2 7,F9.2,T47,F9.2,T67f F9. 2 , T 8 7, F9. 2 , T 10 7 , F9. lift / ) 

WRI TE ( 6,70 

70 FORMAT ( 1H1 , 9 ( / I ,T54, ' RAM PHYSICAL D I MENS I ONS • /T45 , 

1 « { I NCHES EXCEPT WHERE OTHERWISE SPECIFIED)'// > 

WRITE (6, 7 1 )DI A11,DIA12,LN10,LNU,LN12,SPACE1,THC11,THC12,WT1 


■^i 

ON 



71 


FORMAT! T62 » ' MEMBER 1 * /T5 , * 0 1 A22 ' , T20 , 1 01 A1 2 * » T 35 , • LN ID ' , T 50 , 

1 ' L N 1 1 • » T 6 4 , ' L N 1 2 * » T 7 6 , ' SPACE1* f T9 3 t • THCl 1 *, T104, • THC12* , 

2T 1 16, ’WEIGHT ( LB) •/ T4,F6.3 , T 19 , F6.3 ,T34t F 6. 3 ,T49, F 6. 3 , 

3T63, F6*3*T76,F6.3,T89 f F6.3 f Tl33, F6.3,Tii8,F7.3//) 

WRITE (6,72) DI A21»DI A2 2, DI A23 ,LN2 1 , THC21 ,THC22, WT2 

72 FORMAT! T62, •.MEMBER 2 • /T 1 6 , • DI A2 1 * , T32 , * DI A22 * , T48 , • 0 I A23 • , 

1 T64 , * LN21 * ,T79, •THC21* ?T96 * • THC22 * , T 1 1 1 , • WE l GHT ( LB ) • / 
2T15,F6.3,T31,F6.3 ,T47,F6. 3,T63 , F 6. 3 ,T 78, F 6, 3 , T9 5 , F 6. 3 , T 1 13 , F 6. 3/ / ) 

WRITE (6,73 )DI A31 , D I A32 , LN3G , LN31 , LN32 , THC3 1 , THC32 , WT3 

73 FORMAT ( T62 » * MEMBER 3 • /T1 1 , ' 01 A3 , T26, • DI A32 ' , T41 , • LN30 • , 

IT 56, • LN31 • » T7 1 1 ' LN32* ,T86 , • THC 31* T 10 1 , • THC32 • , T 1 13 , « WE IGHT ( LB ) • / 
2T7,F9.3 ,T22,F9.3,T37,F9.3,T52,F9.3,T69,F6.3,T85,F6.3,T97,F9.3, 
3T115,F6.3//) 

WRITE(6,74)DI A41 ,DIA42,LN41 , THC4 1, THC 42 , WT4 

74 FORMAT { T62 ,' MEMBER 4 • / T2 5, • D I A41 » , T40 , • 01 A42* ,T56, • LM41 • , 

1T70* • THC41 1 ,T85, # THC42* , T99 ,• WE IGHT ( LB ) • / 

2T24,F6.3,T39,F6.3,T55,F6.3,T69,F6.3 ,T84,F6.3, TIOi, F6. 3 //) 
WRITE(6,75)DIA52,LM50,LN51,LN52,WT5 

75 FORMAT ( T62 , • MEMBER 5 f /T33, *DI A52' » T49, • LN50 • ? T64, • LN51' , 

1 T78 , «LN52' , T90 , * WE IGHT ( LB ) • / 

2T32 , F6. 3 , T4R , F 6. 3 , 763 ,F6. 3 ,T76 ,F6# 3, T92,F6.3 //) 

WRITE (6,76) WTGAS,WTERR,WTMEM,WE I GHT 

76 FORMAT ( T56, 'WEIGHT SUMMARY! LB ) • / 

1T17, • GAS WEIGHT* ,T4fc, 'WEIGHT-ERROR* , T72, • WE I GHT-M EMBER S' , 

2T102, 'TOTAL WEIGHT* / T 18 , F6. 1 ,T 44, FI 0. 1 , T7 2, F 10. 1 , T 10 1, F 1 G. 1 / // / ) 
WRITE(6,80) 

80 FORMAT( 1H1,18( /) , T54 , ' SPECI FI C GAS CONDIT IONS' 2( /) ) 

WRI TE ( 6 , 8 1 ) MA SSS I * PRESS I ,TEMFLM,vOLS 

81 FORMAT! T56 , ' STAG* CHAM. ,T IME-O. ' / T27, * MASS ( LBM ) ' ,T46, 

1 ' PRES. ( LBF / IN**2 ) ' , T72 , » TEMP. ( F ) • , T9 3 , ' VOLUME ( IN** 3 ) ' / 
2T28,F6.3,T47,F10.1 ,T69,F10.1 t T92,F10. 1////) 

TEMPS(NSTEP$)=TEMPS(NSTEPS)“46C. 

WRITE(6,82)MASSS(NSTEPS) , PRESS F, PRESS ( NSTEPS ) , TEMPS! NSTEPS ), VOLS 

82 FORMAT( T53,'STAG. CHAM. ,TI ME=TFINAL' /T5, ' MASS( LBM ) ' , T20, 

l'MIN. PRES. FOR CHOKED FLOW (LBF/ I N**2) ' , T65, ' ACT. FINAL PRE S ( LBF/ 




2 IN** 2) * rT98, 'TEMP. ( F ) * , T1 1 3 , * V OLUME ( IN**3)*/ 

3T6,F6. 3,T33 r F10. 1,T71,F13.1»T95,F10.1,T112,F1G.1//) 

TEMPS (N STEPS) =TEMPS (NSTEPS )+460. 
TEMPR(NSTEPS)~TEMPR(NSTEPS)-460* 

WRITE <6, 83) MAS SRC NSTEPS) , P RE SR F , TEMPR ( NSTF PS ) t VOL R F 
F 0 R M A T ( T 5 5 t ' R A M CHAM. t TI ME= TF I NAL ' / T28, ' MASSC L BM ) • T46, 
1 'PRES* (LBF/IN**2) ' T72,'TEMP* IF)* ?T92 , 'VOLUME ( IN**3 ) »/ 
2T29*F6. 3#T47, FIO. 1 tT69,F10. 1 ,T91 , F 1 0 * 1/ // / ) 

TEMPR (NSTEPS) =TEMPR(NSTFPS ) 4- 460* 


COMMENT ********** PRINTOUT OPERATIONAL HISTORY ********** 


WRITE!6,90) 

90 FORMAT! 1HI,15!/) ,T59, » HI STORY ( DYNAM I CS ) • // 

1T7,'TIME FROM ZERO ( SEC > ' , T 30 , 1 TI ME BEG. MOT. (SEC)*, 

2T55, • DI SP. ( IN) • » T70 * * VEL. { I N/SEC ) • * T 87 , • FORC E ( BOTT OM , LB F ) • , 

3 T 1 1 0 * 'FORCE ( OUTPUT, LBF) * ) 

DO 91 1=1, NSTEPS 

91 WRITE (6, 92) TIME! I ) ,TlMEAM( I ) ,X< I) , XD! I ) , FRAM! I ) ,FOUT< I ) 

92 FORMAT! T 10 , F 8 . A ,T35 , F8. 4 , T55* F8. A , T70 , F8. 4 , T87, F 1 2. 2 , T 11 0 , F 1 2 . 2 ) 
WRI TE < 6* ICO ) 

100 FORMAT ! 1H1,15!/) »T5 5, 'HISTORY! GAS IN RAM CHAM.)'// 

1 T7, 'TIME FROM ZERO ! SEC) ' *T29 ? ' MASS FLOW RATE! LBM/S EC ) • ,T54, 

2 * MASS ! LBM ) • , To 9, 'PRES. ! LBF/I N**2 ) ' , T 94, • T EMP 1 F ) • , T 1 1 1, 

3 * VOLUME 1 1 N**3 ) * ) 

DO 101 I -1 » NSTEPS 
TEMPR ! I )=TEMPR ! I )-460. 

101 WRI TE 16,1 02 )T I ME! I ) ,MASSDR( I ),MASSR! I ) * PRESR! I ) , TEMPR ( l ) , VDLR ( I ) 

102 FORMAT! T12 , F8. 4 ,T35 , F8. 4, T53, F8. 4 ,T 69, FlO. 1 , T9 1, F 10. 1 , 

1TI09, F10. 1 ) 

WRITE!6,110) 

110 FORMAT! 1H1 , 1 5 ! / ) , T51 , 8 HI STORY ( GAS IN STAG, CHAMBER)'// 

1T5, 'TIME FROM ZERO!SEC)' ,T35,' MASS! LBM ) • , T 53 , • PRES . < LBF/ I N** 2 ) ' , 
2T86, 'TEMP. IF) • f TI06 ,' DENS I TY ! L 8M/I N**2 ) • ) 

DO 111 1=1, NSTEPS 
TEMPS ! I >=TEMPS!I)-46C. 

111 WRITE (6, 112) II ME! I ) ,MASSS! I) , PRESS! I ) , TEMPS! I ) , RHOS! I) 

112 FORMAT! T 1C , F 8 . 4 ,T3 5 , F8. 4 ,T59, F 1 0. 1 , T 83 , F 1 G . 1,T110,F8.4) 

WP I TF ! 6 * 1 20 ) 

120 FORMAT! 1H1 , 1 5 ( / ) , T55 , * HI STORY t G AS IN PORTS)'// 

1T11,«TIME FROM ZERO! SEC) 1 » T7CK *GAS VELOC ITY{ I N/SEC ) • , 

. 2 T 1 0 1 , 'DENSITY! LBM/ I N**3 )' ,T50, ' T EMP! F ) ' ) 

DO 121 1=1, NSTEPS 
TEMPT ( I ) = TEMPT ! I )— 460. 

121 WRITE (6,122) TIME! I ) , TEMPT! I I ,VELGT( I ) , RHOT ! I ) 


122 FORMAT ( T13,F 1 0. A»TA7,F1Q. 1 , T7A, FIG. 1 .T1CA, F 10* 6 ) 

WR I TE ! 6 » 1 30 I 

130 FORMAT ( 1 HI , 1 5 < / I , TA9 , • PORTS AND APPROXIMATE LOCATIONS*// 
l T9 t * PORT NUMBER* ,T29, 'CUMULATIVE ARE A( I N**2 ) * , T59 f 

2 'AREA! IN**2)' ,T86, * DI A! IN I * t T108 , * L OCAT ION! I N > • > 

DO 131 I = 1 » NSTEPM 

131 WRI TE!6,132 INUM! I I » AREACl II t ARE AT ( II , DI AT m » X I I) 

132 FORMAT! TIC t 15 , T33 » FIO . 6 ,T59, F10. 6, T8A, F10. 6, T108 , F 10. 6) 
WRITE ( 6 1 1 AO I 

1 AO FORMAT! 1H1 I 
STOP 



APPENDIX B 


COMPUTER PROGRAM RESULTS 


PYROTECHNIC RAM ANALYSIS 


DIST.BCFT ) 
8:3.300 


TIME PER'C. RAMP 

9C .0 c 


Ki (LBF/ SEC) 

50C0.CC 


K2 (LBF/SEC*** I 
37252.95 


DATA PCI NT-EN E PGY VS WEIGHT CURVE 


ENERGY ( FT-LEF } 
1 3C237.4 


WEIGHT (LB) 
301.5 


CONFIGURATION PARAMETERS 


MLR. BOOSTER ( LBM-FT**2) 
563725152. 


MASS BCCSTER(LEM) 
321000. 


MASS ORB I TER ( LBM ) 
A5CCCC . 


INER. SYS. (SLUG! 
2*79. 


RAM OPERATIONAL CHARACTERISTICS 


TIME OPER.(SEC) 
1. 3C25 


TIME I NCR . (SEC ) 
J. 06 *C 


TIMt KSEC) 
0.0225 


TIME 2( SEC ) 
C.512C 


TIME 3 ( SEC ) 
1.2800 


FURCE 2 ( LBF ) 
2560. CC 


CISP.2I IN) 
0 • 5*134 


VEL.2I IN /SEC ) 
3.17193 


ENERGY 2 (LBF-FT) 
8C.30 


FORCE 3 (LBF) 
1 C C V- CO « CC 


STRCKE(IN) 

29.30183 


VEL . 3( IN/SEC) 
125.8C6C5 


ENERGY 3 (LBF-FT) 
1302C7.** 


FRIC. FORCE (LBF ) 
5CC.C0 


NUMBER OF STEPS ANALYZED 
21 


00 



PYROTECHNIC GAS (TP-Q-3G 27 > PROPER IT I ES 


MOL. WEIGHT 

FLAME TEMP. IF) GAS 

CONST.-R!IN-LBF/LBM-R) 

SP. HEAT RATIO 

CPO( BTU/L9M-R ) 

CVOJBTU/LBM-R) 

21.81000 

1883. 

850.07 

la 1950 

0.2730 

0.2280 

DEN. FACTOR 

PRES. FACTOR 

TEMP. FACTOR 

STAG. PRES. FACTOR 

FLOW FACTOR 

0.62C6 

Q. 5654 

0.9112 

5166983. 

0.4364 


RAM DESIGN PARAMETERS 


DESIGN PRES. -MEM. 2!LBF/IN**2) SAFETY FACTOR DESIGN PRES.-MEMS. 3 AND 4!LBF/IN**2| 

3656.0 1.50 15000. C 


RAM MATERIAL ULTIMATE STRENGTH! LBF/IN**2) YIELD POINT! LBF/IN**2 I DENSITY! LBM/I N**3 ) 

INCONEL 718 18350C. 156C0C. 0.2970 


RAM SIZING AND OPERATIONAL FACTORS 


VOL FAC 

ERROR FACTOR 

THC/PORT DIA. 

INFRTIA FACTOR 

TIMFF/TIMEA 

0.47 

0.20 

20.00 

0.004840 

2.50 


ltifac 

LT2FAC 

LT3FAC 

LT4FAC 

LTSFAC 

D23FAC 

1.10 

1. 10 

1.10 

1.1C 

1. 50 

1.50 



RAM PHYSICAL DIMENSIONS 
IINCHES EXCEPT WHERE OTHERWISE SPECIFIED) 


DIA22 

4.700 


DIA12 LN10 

6. COO 51.953 


MEMBER I 

LNIi LN12 SPACE1 THC1I 

8.000 43.953 1.000 0.650 


THCI2 WEIGHTCLB) 

0.715 175.605 


MEMBER 2 


01 A21 
5.000 


DIA22 
6. 179 


DIA23 
9.2 69 


LN21 

44.953 

THC21 

0.090 


THC22 

0.099 

WEIGHTILBI 

24.854 

DIA31 

2.749 

D I A 32 
3.087 


LN3C 

37.30? 


LN31 

29.302 

MEMBER 3 

LN32 

8.000 


THC31 

0.169 

THC32 
0. 185 

WEIGHTILB1 

17.120 


DIA41 

5.000 


DIA42 

5.613 


LN41 

8.CC0 

MEMBER 4 

THC41 

0.307 


THC42 

0.337 

WEIGHT(LB) 

17.093 




0IA52 

1.000 

MEMBER 5 

LN5C LN51 LN52 

60.953 5.000 3.000 

WEIGHTILB) 

14.218 


GAS WEIGHT 
2.^ 


WEIGHT SUMMARY ( LB ) 

WEIGHT-ERROR WEIGHT-MEMBERS 

49.8 248. ° 


TOTAL WEIGHT 
3C1.5 



SPECIFIC GAS CONDITIONS 


MASS! LB* ) 

1.412 


STAG. CHAM. ,TIME=0. 

MASSILBMI PRES. (LBF/IN**2) TEMP. IF) 

2.851 15000,0 1883.0 


VOLUME! IN**3) 
378.6 


STAG.. 

M!IN. PRES. FOR CHOKED FLOW (LBF/IN**2> 
6465.7 


CHAM. » TI ME =T FINAL 

ACT, FINAL PRESCLBF/IN**2) 
6480. 1 


TEMP. C F J 
1583.1 


MASS (LB Ml 
1.439 


RAM CHAM. »TIME=TF INAL 
PRES. <LBF/IN**2) TEMP.IF) 

3656.0 1968.7 


VOLUME! IN** 3 1 
812.3 


VOLUME! IN* *3 1 
378.6 


* 


HISTORY! DYNAMICS) 


FROM ZERO! SEC) 

TIME BEG. MOT. ! SEC) 

DISP. ( IN) 

VEL. < IN/SEC) 

FORCE! BOTTOM, LBF ) 

FORCE ( OUTPUT , LBF ) 

0.0225 

0.0 

0.0 

0.0 

500.00 

0.0 

0.0865 

0.0640 

c.0011 

0.0496 

820.00 

320.00 

C. 1505 

0.128C 

070085 

C. 1982 

1140.00 

640. CO 

C. 2145 

0.1920 

0.C285 

0.4461 

1460.00 

960.00 

0.2785 

C.256C 

0.0677 

0.7930 

1780.00 

1280.00 

0. 3425 

0.3200 

0.1322 

1.2390 

2100.00 

160C.0C 

0.4065 

0. 384C 

0.2284 

1.7842 

2420.00 

1920.0C 

0. 4705 

0. 4480 

0.3627 

2.4285 

2740.00 

2240.00 

0.5345 

0.512C 

0.5413 

3.1719 

3C60.00 

2 560.00 

0.5985 

0.5760 

0.8556 

4.1895 

4600.62 

4100.62 

0.6625 

0. 6400 

1.2521 

5.7751 

6750.00 

6250. GO 

0.7265 

C.7C40 

1.7738 

8.1390 

9650.62 

9150.62 

0. 7905 

0.7680 

2.4783 

11.5379 

13460.00 

12960.00 

C. 8545 

0.8320 

3.4416 

16.2795 

18350.57 

17850.57 

0.9185 

0.896C 

4.76G9 

22.7275 

24509.93 

24009. 93 

0.9825 

0.9600 

6.5585 

31.3059 

32140.54 

31640.54 

1.0465 

1.0240 

8.9862 

42.5036 

41459.82 

40959.82 

1.1105 

1.0880 

12.2292 

56.8796 

52700.55 

52200.55 

1.1745 

1.1520 

16.5105 

75.0668 

66109. 88 

65609. 88 

1.2385 

1.2160 

22.0966 

97.7772 

81950.56 

81450.56 

1.3025 

1.28C0 

29.3017 

125.8060 

100499.63 

99999.63 


00 

O' 


HI STORY ( GAS IN RAM CHAM. I 


OM ZERO(SEC) 

MASS FLOW RATEUBM/SECI 

MASS f LBM) 

C.0225 

0.0G21 

0.0000 

0.0865 

C. 000 5 

C.0C01 

0.1505 

0.C006 

0.0001 

0.2145 

C.C0C9 

0.0002 

0.2785 

0,0014 

0.0002 

0.3425 

C. 0021 

0.0003 

0.40 65 

0.0031 

0.CC05 

0.4705 

C.C045 

0.0007 

0.5345 

0.0093 

0.0C11 

0.5985 

0.C218 

0.0019 

0.6625 

0.0457 

0.0C39 

0.7265 

C. 089 1 

0.CC78 

0.7905 

0.1672 

0.0153 

0.8545 

0.3048 

0.0292- 

0.9185 

0.5422 

C.0543 

0.9825 

0.9431 

0.0986 

1.0465 

1.6060 

0.1750 

1.1105 

2. 6826 

0.3041 

1.1745 

4.4114 

0.5184 

1.2385 

7.188P 

C-.8688 

1.3C25 

8.9022 

1.4385 


PRES. (LBF/IN**2> 

TEMP ( F I 

VOLUME ( I N*^3 ) 

18.2 

2600.6 

6.9 

29.8 

2598.4 

6.9 

41.5 

2574.9 

7.1 

53.1 

2522.2 

7.7 

64.8 

2446.2 

8.7 

76.4 

2363.5 

10.5 

88.0 

2238.9 

13.2 

99.7 

2229.3 

16.8 

111.3 

2185.0 

21. 8 

167.4 

2645.0 

30.4 

245.6 

2628.8 

41.3 

351.1 

2495.3 

55.6 

489.7 

2363.7 

75.0 

667.6 

2271.7 

101.5 

891.6 

2201.7 

137.7 

1169.2 

2151.4 

187.2 

1508.2 

2114.1 

253.9 

1917.2 

2083.7 

343. C 

2405.0 

2054.5 

460.7 

2981.2 

2019.6 

614.3 

3656.0 

1968.7 

812.3 


00 

-vi 


* 


I* 


HI STORY ( GAS IN STAG. CHAMBER) 


FROM ZEROl SECI 

MASS(LBM) 

PRES. (LBF/IN**2» 

TEMP. IF) 

DENSITY(LBM/IN**2 ) 

Co 0225 

2.8509 

14999.7 

1883.0 

0.0075 

0.0865 

2.8509 

14999 0 5 

1883.0 

0. 0075 

C. 1505 

2.8509 

14999.3 

1883. 0 

0. 0075 

C. 2145 

2. 8508 

14999.0 

1883.0 

0.0075 

C. 2785 

2.8507 

14998.6 

1883.0 

0.0075 

C. 3425 

2.8506 

14997.9 

1882.9 

0.0075 

Co 4065 

2. 8505 

14996.9 

1882.9 

0.0075 

0.4705 

2.8502 

14995.4 

1882.9 

0.0C75 

0. 5345 

2. 8499 

14993.2 

1882.8 

0.0075 

0.5985 

2. 8490 

14987.9 

1882.7 

0. 0075 

0.6625 

2.8471 

14975.7 

1882.4 

0.0075 

Go 7265 

2.8432 

14951.1 

1881. 8 

0.0075 

0.7905 

2.8357 

14904.0 

1880.5 

0.0075 

0.8545 

2.8218 

14816.8 

1878.3 

0.C075 

C. 9185 

2.7967 

14659.4 

1874.2 

0.0074 

0.9B25 

2.7524 

14382.3 

1867.0 

0. 0073 

1.0465 

2.6760 

139C6.4 

1854.2 

0.0071 

1.1105 

2.5468 

13108.2 

1832.0 

0.C067 

1.1745 

2.3326 

11801.6 

1793. 1 

C. 0062 

1.2385 

1.9822 

9715.4 

1722.7 

0.0052 

1.3025 

1.4124 

6480.1 

1583.1 

0.GC37 


OO 

00 






HISTORY ( GAS IN PORTS I 


TIME FROM ZERO(SEC) 

TEMPIF) 

GAS VELOCITY ( IN/S EC I 

DENSITY! LBM/I'N**3 ) 

0*0225 

1674.8 

28934.1 

0.004674 

0.0865 

1674. 8 

28934. 1 

0.004674 

0.15C5 

1674. 8 

23934.0 

0.004674 

0.2145 

1674. 8 

28934.0 

C. 004673 

0.2785 

1674.8 

28933.9 

C. 004673 

0.3425 

1674.8 

28933.8 

0.0C4673 

0.4065 

1674. 8 

28933.7 

0.004673 

0.4705 

1674.7 

28933.4 

0.004673 

0. 5345 

1674.7 

28933.1 

0.004672 

0.5985 

1674.6 

28932.2 

0. 004671L 

0.6625 

1674.3 

28930.3 

0.004667 

0.7265 

1673.7 

28926.5 

0.004661 

0.7905 

1672.6 

28919.0 

0.004649 

0.8545 

1670.6 

28905.1 

0.004626 

0.9185 

1666.9 

28830. C 

0.004585 

0.9825 

1660.3 

28835. 1 

C. 004512 

1.C465 

1648.6 

28756.0 

0.004387 

1.1105 

1628.4 

28617.6 

C. 004175 

1.1745 

1592.9 

28373.5 

0.003824 

1.2385 

1528.8 

27926.8 

0.003249 

1.3C25 

1401.6 

27019.1 

9.002315 


PORTS AND APPROXIMATE LOCATIONS 


PORT NUMBEP 

CUMULATIVE AREA(IN**2) 

AREA ( I N**2 ) 

0 1 A I IN) 

LOCATION! IN) 

l 

0.000016 

O.OOCD16 

C. 004483 

0.0 

0 

0.000016 

0.0 

0.0 

0.001057 

0 

0. COCCI 6 

9.0 

0.0 

C. 008458 

n 

0.000016 

0.0 

0.0 

0.028547 

n 

0. 00C016 

0-0 

0.0 

0. 067668 

n 

0.000016 

0.0 

o.c 

0.132163 

3 

0. OCC016 

o.c 

0.0 

0.22 8378 

2 

0.000034 

C. 000018 

0.004764 

C. 362656 

3 

0. 000069 

0.000035 

0.006707- 

0.541341 

4 

0.CC0161 

0.000092 

0.010826 

C. 855567 

5 

0.000338 

0.000177 

0.015026 

1.252090 

6 

0.000661 

C .000323 

0.020276 

1.773756 

7 

0,001244 

0.000582 

0.027229 

2.478330 

8 

0.002279 

0.0010 36 

0.036313 

3.441604 

9 

0.C04095 

0.001816 

0.048085 

4.760860 

10 

0.007249 

0.003154 

0.063368 

6. 558531 

11 

0.012731 

0.005482 

0.083546 

8.986238 

12 

0.022452 

D. 009721 

0.111254 

12.229156 

13 

0.0 4065 9 

0.018207 

0.152254 

16.510498 

14 

0. 079218 

0.038559 

0.221574 

22.096558 


